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Abstract — The newly emerging theory of compressed sensing 
(CS) enables restoring a sparse signal from inadequate number 
of linear projections. Based on compressed sensing theory, a 
new algorithm of high-resolution range profiling for stepped- 
frequency (SF) radar suffering from missing pulses is proposed. 
The new algorithm recovers target range profile over multiple 
coarse-range-bins, providing a wide range profiling capability. 
MATLAB simulation results are presented to verify the proposed 
method. Furthermore, we use collected data from real SF radar 
to generate extended target high-resolution range (HRR) profile. 
Results are compared with 'stretch' based least square method 
to prove its applicability. 

I. Introduction 

The range resolution of a radar system is determined by the 
bandwidth of the transmitted signal. Stepped-frequency (SF) 
pulse train obtains large signal bandwidth by linearly shifting, 
step-by-step, the center frequencies of a train of pulses. It 
is widely used in high-resolution radar systems and well 
documented in the literature ffl.lZl. In SF radars, the 'stretch' 
processing method 0, based on inverse discrete Fourier 
transform (IDFT) technique, can acquire high-range resolution 
(HRR) profiles with narrow instantaneous bandwidth and 
low system complexity. However, SF radar suffers greatly 
from missing-pulse problems due to interference or jamming 
impinging on the receiver, since SF technique occupies a large 
bandwidth. While some pulses missing are and hence must be 
discarded, the IDFT based stretch processing will inevitably 
leads to high sidelobes, thus undermining the profiling quality. 
Various methods have been proposed to interpolate the missing 
data (see [3| and reference therein). Theoretical analysis and 
experience indicate that the longer the signal interpolation 
length is, the larger the interpolation error is. If the missing 
pulse number becomes bigger, the performance of existed 
method will reduce rapidly 0. 

Besides the missing-pulse problem, SF radar suffers from 
'ghost image' phenomenon. This problem, mainly caused by 
range ambiguity among adjacent 'coarse-range-bins', is deli- 
cately addressed in [4], where the author solved the problem 
by least square (LS) technique. But this method is applicable 
on the assumption that full pulses are well received, and the 
foundation of it is still IDFT. Therefore, missing-pulses also 
deteriorate the profiling results. Recently, the new emerging 
theory of compressed sensing (CS) [5|,[6| that achieves high 
resolution has been widely used in radar applications Q. 
The main advantage of this theory is that, with sub-Nyquist 



samples, sparse signal can still be reconstructed perfectly. CS 
theory was introduced in the signal processing for SF radar 
by Sagar Shah et al. (SJ. With reduced number of transmitted 
pulses in one coherent processing interval (CPI), their method 
provides super-resolution ability in both range and Doppler do- 
main. It also indicted that missing-pulse problem can be solved 
with their method. However, they only discussed profiling 
range of only one coarse-range-bin, limiting their application 
on narrow-range-gate profiling. 

This paper introduces a new profiling algorithm for SF 
radar with missing pulses, and the profiling range gate extents 
for multiple coarse-range-bins. We focus on profiling of a 
stationary object. Unavailable data from missing pulses are 
discarded; sparse recovery is used to obtain extended synthetic 
range profile. We demonstrate that new algorithm can solve the 
missing-pulse problem, it also has a wide profiling range gate. 
The remainder of this paper is organized as follows. In Section 

II, the signal model of HRR profiling for SF radar is stated. 
CS based profiling with missing pulses is described in section 

III. Simulation results are presented in section IV. Section V 
concludes the new approach. 

II. System Model 

In SF radar, a pulse train of N pulses are transmitted with 
stepped carrier frequencies. For the nth pulse, the carrier 
frequency is /„ = f c + nAf, where f c is the initial frequency 
and Af the frequency step. The complex profile of the 
measured scene can be represented by system function H{td), 
as has been derived in is the time domain variable, 

and H(td) describes the complex reflectivity of measured 
scene corresponding to time delay td- For the convenience 
of signal modeling and derivation, it is assumed that one 
target falls in the range gate [Rq,Rq + D] over the whole 
coherent processing interval, where Rq — cQ/2Af, and 
D = cL/2Af (Q,L are nonnegative integers and c is the speed 
of light). In the 'stretch' processing [l j, the range resolution 
is c/2NAf 0. Choosing this resolution as the sampling 
period, the pth high-resolution range cell, which represents 
the complex amplitude of the scatterer located in the range 
R Q + (cp/2NAf), is written by h p = H(p/NAf). Thus, 
the HRR profile of the target can be expressed by the vector 
h = [h ,h 1 ,,h NL ^ 1 ] T . 

A target response matrix (TRM) was used to organize the 
echo signal of the pulse train. The TRM contains N rows and 



S columns. The nth row consists of S uniformly sampled time- 
domain data from the baseband echo signal of the nth pulse 
(If the transmitted baseband waveform is pulse compressing 
waveform, the 'baseband echo signal' refers to the pulse- 
compressed echo signal). The elements in the same column are 
baseband samples of the same coarse range cell. The column 
number is S = 2D/cAt, where At is the sampling interval. 
The TRM of a target can be denoted by 
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Here, E n (r) is the baseband echo signal of the nth pulse, 
and t is the baseband sampling instant. As derived in [1], the 
baseband echo signal from stationary target is 
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where Rx(t) is the baseband pulse shape and u n (r) is 
additive noise. 

In the 'stretch' processing method, the IDFT is applied to 
each TRM column to form a HRR profile in one coarse-range- 
bin fl]. Missing pulse problem means data from some rows of 
TRM are not available. If the missing pulse number is large, 
profiling quality is greatly decreased using IDFT. Our new 
method solve this problem by sparse recovery based on CS 
theory, which can provide a better profiling quality. Based on 
the observation that discrete system function vector h is sparse, 
we propose a new scheme for HRR profiling based on sparse 
recovery in the next section. 

III. HRR PROFILING WITH MISSING PULSES 

We now introduce the new CS based HRR profiling method, 
on condition that some pulses are missing. Suppose only M 
pulses (M < N) from N transmitted carrier frequencies are 
valid, that the carrier frequency of the mth valid pulse is F m — 
fc + C m Af, where m is an integer between and M — 1, C m 
is an integer between and N — 1. Substituting pulse number 
index n in equation (2) by C m , we derive sample output for 
the mth pulse at sampling instance r 



NL — l 



p=0 



c NAf 



(3) 



We rewrite (3) in vector multiplication form: 

Ec m {r) = <f(C m , r)h + u m (r). (4) 

^{Crm t)is a row vector of length N x L, the pth element of 
the vector is 
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Deleting the invalid data in the TRM, 
decreases to M. 



the row number 
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The new TRM includes all available information we received. 
Note that each element of the TRM is a linear projection of 
system function h. By vectorizing this matrix, we may write 
the following equation 

Y = vec(E) = $h + U. (7) 

The observation vector Y is of length M x S. Matrix $ is the 
projection matrix of M x S rows and NxL columns, each row 
of $ is corresponding to an observation. For an instance, the 
row corresponding to pulse number C m and sampling instance 
sAt is ip(C m , sAt). U the noise vector for all observations. 
We have established a linear projection for complex profile 
h. While pulses are missing, M < N holds, and inequality 
MN < SL holds. Therefore, (7) becomes an underdetermined 
equation. According to CS theory, recovering a sparse signal 
from insufficient observation is possible by l\ minimization 
0: 



s.t. \\Y - $h|| 2 < e 
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where h is an reconstruction of h and e is an estimation error 
that is determined by received signal noise. 

IV. Results 

We show some primary results of simulation. The HRR 
range profile of a real aircraft (Fig. 1(a)) was measured by 
a wideband C-band chirp radar. The chirp bandwidth was 
512MHz, providing a range resolution of about 0.3m. This 
measured range profile is used as the scatterer truth. For 
SF radar simulation, 32 LFM pulses are transmitted in a 
coherent pulse train. The frequency step size is 16MHz; and 
the total effective bandwidth is NAf =512MHz. In each 
pulse, single-pulse bandwidth is 24MHz. The sampling rate f s 
equals single-pulse bandwidth. The profile range gate covers 
12 coarse-range-bins. We simulate the missing pulse condition 
by discarding data received from randomly selected 12 pulses, 
the left 20 pulses are valid. White Gaussian noise was added 
to the received data, SNR is approximately 15dB. 

A. Simulated Data 

The results of simulation data obtained via different meth- 
ods are compared in Fig.l. Fig. 1(b) show the result obtained 
from LS method [4], Fig. 1(c) demonstrate the result from new 
approach. From which it can be noted that LS method has 
created high sidelobe, and the result by new method is more 
similar to original target range profile. 

To analyze the profiling results of the two methods quan- 
titatively, we measure the similarity between the simulated 
target and the reconstruction profile by normalized cross 
correlation. Similarity equals 1 means perfect reconstruction. 
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Fig. 1. Comparison of result obtained via different methods for simulative 
data, (a) Model of scatters, (b) LS result of missing data, (c) Sparse recovery 
of missing data. 
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Fig. 3. Comparison of result obtained via different methods for real radar 
data, (a) Full pulse recovery, (b) LS result of missing data, (c) Sparse recovery 
of missing data. 
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Fig. 2. Similarity measure. 



Fig. 2 illustrate the comparison. We increase the number of 
missing pulses from to 20. The line marked by "A" denotes 
the similarity by LS method, and line marked by "□" denotes 
the similarity by sparse recovery. Sparse recovery has an 
obvious advantage over the LS counterpart. 

B. Real Radar Data 

We use real radar data obtained from SF radar. An ex- 
periment was carried out in a wide and fiat field. A single 
metal reflector was placed 1010m away from the radar an- 
tenna. Experimental data of I/Q channels was collected from 
the baseband of the radar receiver. All parameters in the 
experiment were equal to the simulated data. We discard 12 
pulses randomly to simulate the missing data condition. New 
approach is applied to the missing data. Profiling results are 
compared to IDFT based LS method. 

Fig. 3(a) shows the profiling result of LS method with full 
data. Fig. 3(b) and 3(c) compare the profiling results to the 
missing data by LS method and the new method respectively. 
LS method exhibits high sidelobe as predicted, while the 
profiling result by new method is similar to full data profiling. 



Sparse recovery outperforms LS using real radar data. 

V. Conclusion 

The application of sparse recovery in extended HRR profil- 
ing for SF radar is illustrated. The simulated data and real data 
experiments prove that the proposed method is an appropriate 
tool to deal with missing data problem. Profiling quality of the 
new method has an obvious advantage over IDFT based least 
square method, if some pulses are missing. Moreover, it can 
profile multiple coarse-range-bins simultaneously, indicating 
a wide profiling range. The profiling result is not corrupted 
by ghost images. Further work should consider reducing 
computational load for real-time implementations. 
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